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Abstract 

The geometry of rotations in dimensions 3, 4, and 5 is discussed 
using the matrix exponential map. Exphcit closed formulas for the 
exponential of an antisymmetric matrix, as well as the logarithm of a 
rotation, are given for these dimensions. 

1 Introductory remarks 

For representing rotations in three dimensions, one usually considers two 
choices: Euler angles, or quaternions. For aerospace navigation, Euler angles 
are a natural choice, since the familiar maneuvers of yaw, pitch, and roll are 
coordinate axis rotations. However, from a global point of view (as opposed 
to the local viewpoint of a pilot), Euler angles are unnatural and awkward. In 
particular, the decomposition of a rotation about a noncoordinate axis into 
Euler angles is not unique and depends on the sequence of axis rotations. 
On the other hand, quaternions can be used to represent a rotation in a 
more convenient manner: rotation about an arbitrary axis can be written 
in a simple form, and given a rotation, the extraction of the rotation axis 
and angle is a relatively straight-forward process. The price paid for this 
convenience, of course, is the introduction of a less familiar Mathematical 
apparatus. 

In less common use is a third choice: the representation of a rotation by 
a linear map, or matrix. Fundamentaly, the Euler angle representation is a 
matrix representation, but one only uses the representation of coordinate axis 
rotations, which have a very simple form. In some texts, one may encounter 
a matrix formula for rotation about an arbitrary axis, but usually this is 
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presented in the form of a rather unappeahng and unyieldy 3x3 matrix. In 
this article, we will recast this matrix in a more appealing form by using some 
simpler linear maps. As a natural byproduct of this approach, we will obtain 
a simple way of extracting the rotation axis and angle, but using a more 
familiar apparatus than quaternions. Moreover, we will be able to extend 
the constructions to higher dimensions. While one may use quaternions to 
represent rotations in four dimensions (see [5]), one cannot go beyond. 

2 Rotations in three dimensions 

A rotation R in three dimensional space is specified by a (real) scalar, the 
rotation angle 6, and a (real) unit vector u, the rotation axis. The effect of 
R on any vector v is given by Rodrigues ' Rotation Formula: 



Geometrically, if we write v in terms of its components parallel and perpen- 
dicular to u (respectively): v = vy + v_l, then i?V|| = vy, and Rv± is the 
rotation of in the plane orthogonal to u by 6' radians. It is important to 
note that this is counterclockwise rotation about u; indeed, using the stan- 
dard three-dimensional vector identity u x (v x w) = (u ■ w)v — (u • v)w, 
one computes x (u x v) = |v^pu, which is a positive multiple of u. 

2.1 Rotation as a linear map 

We may express equation ([1]) as a 3 x 3 matrix, or equivalently, a linear map. 
Let proj^ denote the linear map proj^{v) = (u ■ v)u (we are still assuming 
that |u| = 1); that is, proj^ is orthogonal projection onto u. In addition, let 
Au denote the linear map Au(v) = u x v. Explicitly as matrices. 



i?v = cos^^ V + (1 — cos0)(u ■ v)u + sin^u x v. 



(1) 




and A 



u 




where u = {ui,U2,U3). Using these maps, equation ([T]) becomes 



R = cos 9 1 + {1 — cos 9) proj^ + sin 9 Au. 



(2) 



Here / denotes the identity map: Jv = v. 



2 



2.2 Finding the angle and axis of a rotation 

Conversely, given a rotation R, we may deduce its rotation angle and axis. 
For the angle, we take the trace of R in equation (|2]) to get tr{R) = 1+2 cos 9. 
Thus 

6 = cos'^{tr{R) -1). (3) 

For the axis, we make use of the fact that / and proj^ are symmetric, while 
Au is antisymmetric, so that equation (j2]) also yields 

= sin^Au, (4) 

where i?* denotes the tranpose of R. The components of u can be read 
directly off of A■^^ = {R — R*)/ 2 said (or more practically, since 2sin6'Au = 
A2sineu; wc cau just read off the components of 2 sin 6' u from R — i?*, and 
then normalize to get u). 

Two remarks are in order. First, the rotation angle 9 determined by 
equation ([3]) is in the range < ^ < tt, whereas there is no such restriction 
in equation ([2]). If the original rotation angle is greater than vr (but less 
that 27r), then the axis determined from (jl]) will be opposite of the original 
rotation axis. Second, we can only use equation (j4]) to determine the rotation 
axis if < ^ < TT. In the case 6' = 0, the rotation is the identity. And in the 
case 6 = 7T, the rotation is reflection in the plane orthogonal to the vector u, 
which we can deduce up to sign using proj^^ = |(/ + i?). 

2.3 Rotation as an exponential 

Recall that the Taylor series expansion of the exponential function is given 
by = ^fc>o3^'^/^'5 which converges for all real numbers x. Similary, the 
exponential of a square matrix A is defined by the formula 

oo ^ 

exp(A) = 5^-A^ 

k=0 

Here we set A° = /. It can be shown that this summation converges for any 
matrix square A (see for example, |T]). 

It is known from the theory of Lie groups that every rotation is the ex- 
ponential of an antisymmetric matrix. We can verify this directly in three 
dimensions. Let u be a three-dimensional unit vector. Starting from the 
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matrix representation of Au, or from vector identities, one computes that 
= —(J — proj^) and A^ = — Au. In particular, — A^ is orthogonal projec- 
tion onto the plane orthogonal to u. It follows that A^^ = (—1)^(1 — proj^) 
for all p > 1, and A^^"^^ = (— l)^Au for all p > 0. We then compute 

k=0 p=l ^ ^' p=0 ^ ' 



1 + 



02p 



[I - projj + 



_p- (^py- 

I + {cos 9 — — proj^) + sin^^Au 



E 

,p=0 



{2p + l)\ 



02p, 



by virtue of the Taylor series expansions for sine and cosine. This last ex- 
pression is just a rearrangement of equation (|2]). That is, exp(6'Au) is coun- 
terclockwise rotation by 6 radians about u. 



3 Simple rotations 

We now extend the above constructions used in three dimensions to any 
dimension. To do this, we define two special linear maps. If u and v are 
two vectors in n-dimensional space, their outer product is the linear map 
u ® V such that 

(u v)(w) = (v ■ w)v 

for all vectors w. This is a generalization of the projection operation: if 
|u| = 1, then proj^^ = u ® u. The wedge product of u and v is defined as 

uAv = u(g)v — v(g>u. 

The wedge product will serve as our generalization of Au; in fact in three 
dimensions, we have u A v = — Auxv 

3.1 Fundamental computations 

The following properties of the outer and wedge products are readily verified 
from the definitions. 

Lemma 3.1. For all vectors a, b, u, v, we have 
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i. tr{u ® v) = u • V 

ii. (a (g) b) (u (g) v) = (b ■ u) (a (g) v) . 

in. |u A vpj)roj^y = |vpu ® u — (u ■ v)(u (g) v + v ® u) + |upv ® v. 

iv. (u A v)^ = — |u A vpjjrojuv 

w. (u A v)^ = — |u A vpu A V 

where |uAv| = a/Iu^IvP — (u ■ v)^, andproj^^ is orthogonal projection onto 
the plane spanned by u and v. □ 

We remark that, as the notation indicates, |u A v| is the matrix norm 
of u A V. For antisymmetric (or symmetric) matrices A and B, one usually 
defines their inner product a.s A ■ B = ^tr{AB^). In this sense, |u A vp = 
(u A v) ■ (u A v). In particular, u A v = if and only if |u A v| = 0. Note 
that we may also write |u A v| = |u| |v| sin0, where is the angle between 
u and V. 

Using lemma I3.1[!ivl) and (jvj) , we may again sum the Taylor series as we 
did in section [231 to compute the exponential of u A v. 

Theorem 3.2. For all vectors u, v with |u A v| = 1, and any real scalar 9, 
exp(6' u A v) = / — (1 — cos 0) proj^^ + sin 6* u A v. □ 

3.2 Rotation in a two— plane 

In any dimension greater than one, a (linear) two-plane is spanned by two 
linearly independent vectors. We will call a linear map corresponding to 
rotation in a two-plane, and that fixes the orthogonal complement, a simple 
rotation. 

Lemma 3.3. Given vectors u, v, if a = au + /3v and b = 7U + (5v for 

scalars a, (3, 7, 6, then a A b = {a6 — fi'-^) u A v. 

Proof. Use the linearity and antisymmetry of the wedge product. □ 

Lemma 3.4. The vectors u, v span a two-plane if and only ^/u A v 7^ 0. 
Moreover, the vectors a, b span the same two-plane if and only if a Ah is a 
nonzero multiple o/u A v. 
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Proof. The first statement follows from the remarks after lemma 13.11 The 
necessity of the second follows from the previous lemma; for the sufficiency, 
use lemma [3?Tlliv|) to show that proj^^y, = proj^^. □ 

An orientation for the two-plane spanned by u, v is a choice of ordering 
of the vectors, say (u, v). By lemma [231 a pair of vectors (a, b) from this 
two-plane has the same orientation if a A b is a positive multiple of u A v. 

Theorem 3.5. // |u A v| = 1, then exp(^u A v) is simple counterclockwise 
rotation by 9 radians in the oriented two-plane spanned by (u, v). 

Proof. By lemma 13.41 we may assume u and v are mutually orthonormal. 
Using theorem 13. 2[ we then compute exp(^u A v)u = cos6'u — sin^ v, and 
exp(^ u A v) V = sin 6u + cos 6 v. Moreover, if w is orthogonal to both u and 
V, then exp(6' u A v)w = 0. □ 

3.3 Logarithm of a simple rotation 

From theorem 13.51 we see that every simple rotation R can be expressed as 
the exponential of the wedge product of two vectors, say R = exp(u A v). In 
this case, we refer to u A w as a logarithm of R. Logarithms of rotations are 
not unique, as the sine and cosine functions are periodic. However, for simple 
rotations by angles less that tt, there is a canonical choice for logarithm. 

Theorem 3.6. If R is a simple rotation in dimension n, then the angle of 
rotation is given by 9 = cos~^ |(tr(i?) —n + 2). Moreover, if0<9<7f, then 
R = exp{f), where f = {e/2sme){R - R'). 

Proof. From lemma 13. tr{proj^^) = 2 and tr{u A v) = 0. Thus from 
theorem 13.21 tr(exp(^^u A v)) = n — 2(1 — cos 9). □ 

If < ^ < TT, we write log(i?) = /. If ^ = 0, we simply take log(-R) = /. 
And if ^ = vr, then we are only able to deduce the projection onto the two- 
plane: proj^^ = ^{I — R); R is multiplication by —1 in this two-plane. 

4 Decomposing an antisymmetric matrix 

In two and three dimensions, every rotation is simple. In higher dimensions, 
one must compose simple rotations to obtain a general rotation. If the two- 
planes of the constituent simple rotations are mutually orthogonal (from Lie 
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Theory, this can always be done), we may obtain an exphcit formula for the 
composite rotation. We will do this in the next section, although we outline 
the general procedure here. 

Lemma 4.1. The two-planes spanned by the pairs u, v and a, b are orthog- 
onal if and only z/ (u A v)(a A b) = 0. 

Proof. Since u-a = u- b = v- a = v- b = 0, lemma I3.1l1ii|) implies that 
(uAv)(aAb) = 0. Conversely, if (uAv)(aAb) = 0, then item lemma [3miiv|) 
implies that proj ^^proj = 0, and the two-planes are orthogonal. □ 

Suppose that fj, I < j < m, are wedge products with mutually orthgonal 
two-planes, lemma WTLl then implies that (^^- fj)'^ = Y2j fj fo^' k > 0. By 
summing the Taylor series, we thus find exp(^^. fj) = I + J^ji^'^vifj) ^ -^)- 
Theorem 13.21 can be then used to obtain an explicit formula for the latter. 

If we are able to write an antisymmetric matix / as a the sum of wedge 
products fj whose two-planes are mutually orthogonal, we will call the de- 
composition / = Ylj fj orthogonal decomposition of /. By dimen- 
sional considerations, there can be at most \n/2\ (where \_x\ denotes the 
integer part of x) nozero summands in the orthogonal decomposition of a 
n X n antisymmetric matrix. 

4.1 Existence of orthogonal decomposition 

Up to this point, we have used only elementary considerations; however the 
the following lemma makes use of the spectral theorem for real nxn antisym- 
metric matrices A: there exists an orthogonal matrix P such that A = PBP^, 
where B has 2x2 blocks of the form ( o ) along the diagonal; if n is odd, 
B has at least one additional row and column of zeros. See [2]. 

Theorem 4.2. Every real antisymmetric matrix, is the sum of wedge prod- 
ucts whose two-planes are mutually orthogonal. 

Proof. Let B and Q be as above; and let e^ denote the n-dimensional vector 
whose j-th component is unity, and all other components are zero. We 
may write B = J2j ^j^j-i ^ ^j+i- Thus B is evidently the sum of wedge 
products with mutually orthogonal two-planes; and since Q is orthogonal, 
A = QBQ^ = XjiQej^i) A (Qej+i) is also. □ 

Corollary 4.3. Every real antisymmetric nxn matrix, with n = 4,5, is the 
sum of at most two wedge products whose two-planes are orthogonal. □ 
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4.2 Finding an orthogonal decomposition 

Theorem 4.4. Let f be a real n x n antisymmetric matrix, with n = 4,5, 
then 

A = tr{f) - \tr\f) and el = -\tr{f) 

are nonnegative. Moreover if A 0, then we have the orthogonal decompo- 
sition f = f+ + f- where 

are wedge products with f+f- = and |/±| = 6'± = \/6±. This decomposition 
is unique up to reordering of summands. 

Proof. By corollary 14. 3[ we may write / as the sum of wedge products / = 
f^ + f-, with /+/_ = 0. However from lemma [3?T| = —\f±\'^proj^, where 
proj_^ denotes orthogonal projection onto the two-plane of the wedge product 
f±. Therefore, we have = —\f+\'^proj^ — \f-\^proj and by taking the 
trace, we get the equation (*) + = — | tr(/^). Similarly, as the 
square of a projection is itself, one computes the trace of f^ to obtain the 
equation: (**) |/+|'' + |/_|^ = | tr{f'^). Settings = and ?/ = these 
two equations give us a system of two algebraic equations in the variables 
X, y. Using (*) to eliminate y in (**), we arive at the quadratic equation 

+ ^tr{f'^)x + ^tr^[f'^) — \tr{f^) = with discriminant A, and whose 
solutions are x = 61- With possible reindexing, we may take x = and 
y = el. Note that (*) and (**) also imply that A = 2(|/+|H - + 
= - > 0. In addition, = > 0. 

To obtain the formulae for f±, we use lemma [XTTrvj ) to get f^ = f\ + f^ = 
— — elf^. This leads to the Vandermonde system 




The determinant of the 2x2 matrix on the left is el—e^ = —\/A. Since we 
are assuming A > 0, we may invert the system to get the stated (necessarily 
unique) solutions for /+,/_. □ 

In the case when A = (so that 6*^ = 0^ ), an orthogonal decomposition 
exists, but will not be unique. This stems from the fact that an eigenspace 
of dimension greater than unity has no natural choice of basis. 
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5 Rotations in four and five dimensions 



As we mentioned previously, every rotation is necessarily the exponential of 
an antisymmetric matrix; see [1]. We use this fact, coupled with orthogonal 
decomposition, to write down a parametrization for all rotations in four and 
five dimensions, and to find logarithms for such rotations. 

5.1 Rotations 



Theorem 5.1. Let f, A, 9± be as in theorem 4-4' U ^ ^ o,'^'^ ^± 7^ 0? 
then 

exp(/) = I+^[Af + Bf + Cf + D/^), 
V A 

where 

A = -L_ _ — 1^ 5- +^ ' ^' 



e_ 9+ ' 9i 9: 



2 02 



. sin 9^ sin 9^ 1 — cos 9^ 1 — cos 9^ 

^ ~— 7. 7. 1 ~ 



9_ 9+ ' 91 91 

sixi X cos 

Moreover if A = and 9^ ^ 0, then exp(/) = / H — / H p . 

9+ 9_^ 

Proof. By theorem iJ we have /+/_ = 0, so that {9lf + f^){9lf + f^) = 0; 
from which it follows that = -9l9lf^ - {9l + 9l)f\ If A ^ and 9+ ^ 0, 
then lemmaOJivD implies proj^ = -{l/9lA){9lf + f^y = {l/9lA){9lf^ + 
/^), and similarly for proj_. The first stated formula for exp(/) then follows 
from theorem 13.21 and the discussion at the beginning of section HJ On the 
other hand, if A = and 6*+ 7^ 0, then 6*+ = 6*-; so that |/+| = |/_|. Thus 
f = fl + fl = -\f+\V+ + /-) by lemma OJvj); whence = -9lf. 
We may then sum the Taylor series for exp(/) to obtain the second stated 
formula. □ 

The case A 7^ and 9+ = cannot occur, since —tr{f^) = |/p > 0; if 
9_ = 0, then / is already a wedge product, and theorem 13.21 applies. There 
are no other nontrivial cases. We remark that in the case when n = 4, A 7^ 0, 
and 9± 7^ 0, we can use the identity projj^ + proj _ = I to eliminate in the 
formula for exp(/). 
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5.2 Logarithm of a rotation 

Theorem 5.2. Let R be a n x n rotation matrix, where n = 4,5. Define 
g = \tr{R^) - \tr'^{R) + \{n - A)tr{R) - \n{n - 6) 

y± = \{tr{R)-n + A)±\^6. 

Then 5 > and — 1 < ?/± < 1. If 6 and y± ^ ±1, then exp(/) = R, 
where f is antisymmetric and has the orthogonal decomposition f = + 
with 

f± = . J yAR - R') -\{R'- R'')), 

2 sm o'-tv 

where 6± = cos~^y±. If 6 = and y+ ^ ±1, then f = ^——{R — i?*). 

2 sin 9+ 

Proof. We know that R = exp(6'+w+ + 6'_w_), where w± are wedge products 
of unit length and whose two-planes are orthogonal, and where < ^± < tt. 
By theorem 13.21 and the discussion at the beginning of section HI 

i? = / — (1 — COS 9^)proj_^_ + sin6'+w+ — (1 — cos9_)proj _ + sin6'_w_, (5) 

where proj_^ = — w^.. Using lemma [XT]fliv|) and (jvj), orthogonality, and the 
double angle formulas sin 26' = 2 sin cos ^ and cos 26' = cos^ ^ — sin^ ^, one 
computes 

R^ = I — (1 — cos 29^)projj^ + sin 29^w^ 

— (1 — cos29J)proj _ + sin20_u;_. 

Set y = cos 6'+ and z = cos 9_. Since cos 29 = 2 cos^ 6*— 1, we get the equations 
(*) tr{R) = 2y + 2z+n — 4: and (**) tr{R^) = Ay'^+Az'^+n — S; and from these 
we get y'^-l[tr{R)-n+A]y+l[tr^{R)-tr{R'^)-2{n-A)tr{R)+n^-7n+8] = 0, 
whose discriminent is S, and solutions are y = y±. Moreover, one computes 
using (*) and (**) that 6 = (?/+ — ?/_)^ and we may take y± = cos9±. Thus 
5 > and -1 < ?/± < 1. 

Let /+ = 9^w^ and /_ = 9-W-. First assume that 6^0 and < 
9± < TT (that is, y± ^ ±1). Equation ([5]) then implies that \{R — i?*) = 
sin6'+/+/6'+ + sin9^f^/9^, and equation ([6]) implies that \{R^ — -R^*) = 
sin26'_|_/+/0_i_ + sin2^^_/_/^^_. Using the double angle formula for sine, we 
may recast these equations as the matrix equation 

1 1 \ /sin^^+/+/^^+\ ^ / {R~R')/2 
2y+ 2y]Uin9_f_/9_ UR^-R^')/2 
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The determinant of the 2x2 matrix on the left is —2\f5 ^ 0, so we can invert 
to get the first solution stated. Second, if 5 = 0, then ?/+ = y_. Therefore, 
equation ([5]) becomes R = I — {1 — y+){proj_^_ + proj_) + sinO^f /9+, which 
implies the second solution, provided that < 6*+ < vr. □ 

In the case when 6^0 and ?/+ = 1 or = 1, i? is a simple rotation, 
and is handled by theorem 13.61 Note that the case 5 7^ and y+ = —1 
is not possible (otherwise ?/_ < —1). Moreover, from equation (E]), we see 
that if y± = —1 (so that 6 = 0), then R is multiplication by —1 in the 
four-plane with orthogonal projection map (/ — R)/2 — although this is 
only meaningful in the case n = 5. The only other nontrivial case occurs 
when 6 and y- = —1; in which case R = exp(/+) — 2proj_, where 
/+ = (^^+/2sin^+)(i? - i?*) and proj_ = (exp(/+) - R)/2. 

6 Concluding remarks 

The formulas given for rotations and their logarithms in three dimensions 
seem to be well-known, although surprisingly obscure — most likely due to 
the predominant use of Euler angles. On the other hand, I have not seen the 
formulas given here for four and five dimensions previously in print, though 
it strikes me as unlikey that they are original. Quaternionic formulas for the 
two-planes and rotation angles for a four-dimensional rotation, however, can 
be found in [3]. 

The constructions given for the orthogonal decomposition, exponential 
map, and logarithm map may be extended to higher dimensions. However 
for the two-plane rotation angles in dimension n, one must solve a polynomial 
of degree \n/2\, which cannot be done analytically for n > 10. Moreover 
to find the orthogonal decomposition, one needs to invert a \n/2\ x \n/2\ 
matrix. 
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